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A METHOD OF ENHANCING SPECTRAL DATA 

Technical Field 

This invention relates to a method of enhancing 
spectral data such as, for example, optical spectral data 
and mass spectral data, obtained from a spectrometer. 

Background to the invention 

Spectral data comprises a series of peaks and troughs 
which correspond to species or elements present within a 
sample (often graphically represented as a graph of 
intensity versus wavelength, frequency, energy or mass) . For 
the case of optical emission spectra, a sample can be 
excited using various known techniques. The excitation 
causes the energy of atoms to be elevated to a higher energy 
level . As the atoms in the excited sample relax or decay to 
a lower energy level of excitation, photons are emitted 
having a discrete wavelength, thereby producing a series of 
sc-called spectral lines, each line corresponding to an 
energy transition. The energy, and hence wavelength, of the 
emitted photon is dependent on the energy gap between the 
excited and relaxed state of the atom, amongst other things. 
The energy levels arid the gap between excited and relaxed 
states are dependent on the atomic element being excited. 
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Thus, it is possible to deduce the constituent elements in a 
sample by looking at the wavelengths of optical emissions 
from an excited sample. 

A typical arrangement for a spectrometer of this kind 
5 is shown in figure 1, which shows, highly schematically, a 
spectrometer . 10 such as the ARL QUANTRIS ™ spectrometer 
manufactured by Thermo Electron Corporation. Here, an 
excited sample S emits radiation 12 which comprises many 
spectral lines. Optical objectives 14, 16 and 18 

10 respectively, each sample a. portion of the radiation into 
the spectrometer. The radiation passes through entrance 
slits 20, 22, and 24 respectively. The sampled radiation 
then impinges on wavelength dispersing elements 26, 28 and 
30 respectively. In this arrangement, the elements are 

15 reflection gratings, known in the spectrometer art although 
other types of dispersing elements can be used. Each grating 
determines the wavelength of radiation which is reflected 
onto detector arrays 32, 34 and 36 respectively, disposed in 
the focal plane of the grating. 26, 28 and 30 respectively. 

20 Modern optical spectrometers such as the 

spectrometer 10 of Figure 1 use solid-state detectors, such 
as CCD (Charge Coupled Device) or CID (Charge Injected 
Device) which comprise at least one array of photo-detectors 
arranged downstream of a wavelength dispersing element and 
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radiation source containing the excited sample. The array is 
mounted in the focal plane of the dispersing element. 
Because each detector has a finite physical width, each 
detector detects a band of wavelengths which is dependent on 

5 the width of each detector element, the dispersive power of 
the dispersing element, and the distance between the 
detector and dispersing element, among other factors. As a 
result, the resolving power of the spectrometer is limited 
by the number of detectors in the array and by the bandwidth 

0 of each detector. 

Referring to figure 2, a typical spectrum 50 recorded 
with the ARL QUANTRIS ™ spectrometer is shown. The spectrum 
is a graphical representation of the recorded radiation . 
intensity 52 (in arbitrary units) as a function of 

5 wavelength 54 (measured in nanometers) . The sample used to 
create this spectrum is composed of pure iron. As can be 
seen, the spectrum is very complex with more than 6000 
spectral lines being visible on a 8640 pixel CCD array. A 
spectrum from a multi-element sample can be more complex 

0 still, depending on the concentration levels of the elements 
that compose the sample. 

Figure 3 shows a portion of the spectrum of figure 2, 
which is presented as a bar graph 60 of the signal detected 
by each individual detector element . The graph clearly shows 
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that the resolution of the detected wavelengths is limited 
by the size or resolving power of the detector array. Each 
detector element effectively has a width over, which it 
integrates the radiation incident upon it, creating a small 
5 wavelength pass band. The size of this pass band (usually 
quoted in picometers wavelength) depends, amongst other 
things, on the physical size of the pixel detector. 

Several line types appear in figure 3: a single line 62 
in the middle, with a FWHM (Full Width at Half Maximum) of 

10 roughly 2-3 pixels or detector elements; two overlapped 
lines, 64 and 66 respectively, to the left; and an 
unresolved group of peaks 68 to the right. The single line 
62 is not centered on a pixel so that it is difficult to 
determine accurately where the centre of the line is . Not 

15 knowing the central position of the line, together with 

having too few digital measurements across the line, makes 
calculating the intensity of the line difficult, and hence 
deducing the quantity of the element that produced the line 
is also difficult. 

20 Present methods which might increase the peak 

positioning accuracy, in most of the cases, use fitting 
techniques. Fitting techniques, such as Gaussian, Lorentzian 
and polynomial (parabolic) fits have proved unsatisfactory 
mainly because these techniques are not able to provide 
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sufficient accuracy of the peak wavelength of a spectral 
line. Inaccuracy on the peak wavelength is typically 
observed when the peak shape is not ideal, for instance 
asymmetrical or due to overlap with one or several other 
5 spectral peaks. Also, such techniques do not allow accurate 
measurement of the peak intensity, and hence element 
. concentration in the sample. Using such fitting techniques 
to calculate spectral line characteristics (such as line 
maximum position, maximum intensity and peak width, for 

10 instance) requires the raw data to ideally have a perfect 
shape, that is the line should be symmetrical, free from 
interference (which might be caused by proximate spectral 
lines causing overlap) , and have a profile corresponding to 
the fitting curve (Gaussian profile if a Gaussian fit is 

15 being used) for an accurate fit to be performed. It is 

highly unlikely that these conditions will ever be realised 
in real spectral data, for instance distortions in line 
shape may occur due to optical or instrument aberrations, 
spectral line overlap, doublet interferences for example. As 

20 a result, improvement of the spectral line shape is often 
not satisfactory using these known techniques. 

The problems outlined above limit the performance of 
optical emission spectrometers, as well as other kinds of 
spectrometers, for a given cost. 
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Summary of the Invention 

It is desirable to increase the digital resolution in 
an attempt to resolve the problems discussed above. 
5 Additional improvements to spectra are also desirable, such 
as improvements to the signal to noise ratio, and signal 
interpolation. 

Against this background, the present invention provides 
a method of enhancing spectral data, said data comprising M 
10 discrete intensity values within a range of wavelength, 
frequency or mass values, 

said method comprising: 

a) applying a first function to the spectral data to 
obtain an inverse transform of the spectrum, 

15 b) zero-filling said inverse transform, and 

c) applying a second function to the zero-filled 
inverse transform to obtain a spectrum comprising N discrete 
intensity values within said range of wavelength, frequency 
or mass values, and wherein N>M. 

20 In embodiments of the present invention, a spectrum is 

measured in the wavelength, frequency or mass domain (or any 
other related domain such as but not limited to energy) , and 
an inverse Fourier Transform (for example) is applied to the 
data to give a spectrum in the inverse transform domain. In 
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the case of spectra which comprise a plot of intensity 
against wavelength or frequency, this inverse transform 
domain is a pseudo-time domain. Zero-filling and, 
optionally, apodization, and then a Fourier Transform is 
5 applied to this pseudo-time domain data to obtain an 
enhanced spectrum in the inverse frequency (wavelength) 
domain. In the case of a mass spectrum, the inverse 
transform domain is not analogous to the time domain but the 
technique can nevertheless be applied equally to it. In 
10 other words, following Zero-filling and, optionally, 

apodization, the mass spectrum can be reconstituted from the 
thus modified data in the inverse transform domain, by' 
applying a second function such as a Fourier Transform' to 
it. 

15 The overall resolution of a spectrometer is a 

combination- of the digital and the spectral resolution 
thereof. Here, the term "digital resolution" is employed to 
describe the resolution of the signal limited by the 
wavelength or frequency interval, or the mass interval, 

20 between two- discrete consecutive values. In a raw spectrum, 
the digital resolution is thus limited by the bandwidth of 
the pixel and (for some detectors) the dead space between 
pixels. The term "spectral resolution", by contrast, 
describes the optical or mass resolution limits of the 
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optical or ion optical components prior to the detector, 
which may include an entrance slit and the dispersive 
element, for example. These two resolution limits are 
combined when the spectrum is measured and the combination 
5 of the two results in a resolution lower than each 
individual resolution. 

The present invention addresses the digital resolution 
of the spectrometer (the spectral resolution being 
determined by the arrangement and components of the 

10 spectrometer itself) . By manipulating the spectral data in 
the time domain, rather than seeking to interpolate the 
"raw" spectrum (that is, the data in the 
frequency/wavelength/mass/energy etc domain) , several 
advantages accrue. For example, the location of the peaks 

15 (both in terms of intensity and in terms of wavelength or 
other spatial position) can be determined more accurately. 
Integration limits can be set with far greater precision. 
Drift compensation (that is, the drift of the spectrometer 
with time as a result of temperature changes etc) can also 

20 be applied more precisely. 

The first function can be a Fourier Transform function 
which produces an inverse Fourier Transform of the spectral 
data. When the spectral data is a wavelength spectrum then 
the inverse transform is of a time-domain interf erogram 
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type. In other words, the spectrum can be transformed into a 
time -domain- like acquisition by inverse Fourier Transform or 
any transform producing a comparable effect. 

Preferably, the first function is an inverse Fourier 
5 Transform (IFT) . Again when the acquired spectral data 
contains intensity as a function of wavelength, applying 
such an IFT to the spectrum transforms it into a time-domain 
like acquisition, hereafter referred to as pseudo-time 
domain signal or interf erogram. This interf erogram is 

10 somewhat analogous to acquired signals from known Fourier 
Transform (FT) instruments (for instance Fourier Transform- 
Nuclear Magnetic Resonance (FT-NMR) , Fourier Transform Ion 
Cyclotron Resonance Mass Spectrometry (FT- ICR) , Fourier 
Transform Mass Spectrometry (FT MS) , Fourier Transf orm- 

15 Infrared (FT-IR) , and so forth). The skilled person will of 
course understand that, in the case of FT -MS and the like, 
the signals are acquired directly in the (true) time domain 
and are then transformed using an FFT into the frequency 
domain and from there (usually) into a mass or other spectrum 

20 - no inverse transform takes place. 

The second transform stage is a transform function, the 
reciprocal of the first function, which transforms the 
signal back to the spectral representation of the signal . 
Other functions (and their reciprocal functions) might also 
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be used to produce similar transformations (for instance z- 
transform, Hadamard transform) . 

Preferably, the invention further comprises the step of 
apodizing said zero- filled inverse transform prior to 
5 applying the second function. The second function can be 
applied to the apodized zero-filled inverse transform. 
Apodization can be used to improve signal-to-noise ratios of 
the enhanced data. 

Furthermore, when the inverse transform is zero-filled 
10 by a factor Z, N is Z times greater than M. Preferably, Z 
should be in the range of 2 to 10. When Z is greater than 
10, there is a burden on computing the enhanced data. Of 
course, as computational methods advance, values of Z>10 may 
be used to great effect. The ceiling value for Z of 10 is 
15 not considered limiting, and higher values might be used 
without leaving the scope of the invention. 

The present invention also provides a computer program, 
which when run on a computer, carries out the method steps 
described above. 

20 Furthermore, the present invention provides a processor 

configured (a) to receive spectral data from a spectrometer, 
the spectral data comprising M discrete intensity values 
within one of a range of wavelength values, a range of 
frequency values and a range bf mass values; (b) to apply a 
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first function to the spectral data to obtain an inverse 
transform of the spectrum, (c) to zero-fill said inverse 
transform, and (d) to apply a second function to the zero- 
filled inverse transform to obtain a spectrum comprising N 
5 discrete intensity values within said one of said ranges of 
wavelength,, frequency and mass values, and wherein N>M. 

The invention may also extend to a spectrometer 
arranged to generate an array of spectral data comprising M 
discrete intensity values within one of a range of 
10 wavelength values, a range of frequency values and a range 
of mass values, the spectrometer including such a processor. 

- In summary, embodiments of the present invention 
provide some or all of the following advantages: 

(A) Peak position accuracy and precision is improved 
15 allowing correct identification of spectral lines and/or 

accurate and precise calibration of the spectrometer; 

(B) The ability to reveal detailed features of the 
optical spectrum obscured by the limited digital resolution; 

(C) Improved accuracy and precision of quantified 
20 spectral features such as peak height and/or peak area 

(setting the integration limits) ; 

(D) More successful drift compensation, drift being 
caused by optical components shifting due to temperature 
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changes, and/or the source position shifting due to argon 
gas pressure changes within the source chamber; 

(E) Increased speed of analysis of sequential 
spectrometer instruments; 

(F) Reduced array size for the detector, reducing 
cost ; 

(G) Ease of applying methods of improving the signal 
to noise ratio; and 

(F) Ability to apply the method retrospectively to 

spectra taken before the ideas set out herein were 
developed, or before the most recent large array detectors 
were available 

Description of the drawings 

An embodiment of the present invention is now 
described, by way of example, with reference to the 
following drawings, in which: 

Figure 1 shows a schematic diagram of a known optical 
emission spectrometer; 

Figure 2 shows a portion of an optical emission line 
spectrum for Fe; 

Figure 3 shows a portion of the spectrum of figure 2; 

Figure 4 shows the data from figure 3 after shifted 
IFFT function has been imposed on the data; 
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Figure. 5 shows raw data obtained from an optical 
emission spectrometer; 

Figure 6 shows the data of figure 5 once it has been 
enhanced according to a first embodiment of the present 
5 invention; . 

Figure 7 shows the data of figure 5 once it has been 
enhanced according to another embodiment of the present 
invention; 

Figure. 8 shows the data of figure 5 once it has been 
10 enhanced according to a further embodiment of the present 
invention; 

Figure 9 shows another data set representing raw 
spectral data; 

Figures 10 and 11 show the data set of figure 9 after 
15 linear interpolation has been performed; 

Figures 12 and 13 show the data set of figure 9 after 
functions according to embodiments of the invention have 
been performed; 

Figure 14 is a plot of the data of figures 11 and 13 
20 superimposed on the same graph; 

Figure 15 is a plot showing various different 
apodization functions which can be applied to data; 

Figure 16 is a plot of raw data obtained from a pure 
aluminium sample; 
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Figure 17 shows the data of figure 16 after a method 
embodying the present invention has been applied without 
apodization; 

Figure 18 shows the data of figure 16 after a method 
5 embodying the present invention has been applied with, 
apodization; 

Figure 19a shows a raw mass spectrum, with the 
intensity and mass axes both shown on a linear scale; 

Figure 19b shows the raw mass spectrum of Figure 19a, 
10 with the intensity axis on a logarithmic scale but with the 
mass axis on a linear scale; 

Figure 20a shows the data of figures 19a and 19b after 
a method embodying the present invention has been applied 
without apodization, and with the intensity and mass axes 
15 both shown on a linear scale; 

Figure 20b shows the mass spectrum of Figure 20a, with 
the intensity axis on a logarithmic scale but with the mass 
axis on a linear scale; 

Figure 21a shows the data of figures 19a and 19b after 
2 0 a method embodying the present invention has been applied 
with apodization and with the intensity and mass axes both 
shown on a linear scale; and 
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Figure 21b shows the mass spectrum of Figure* 21a, with 
the intensity axis on a logarithmic scale but with the mass 
axis on a linear scale. 



5 Detailed description of embodiments of the present invention 
An embodiment of the present invention comprises a 
method of manipulating digitised spectral data which can 
produce a resultant spectrum that more accurately resembles 
the physical spectrum emitted from the sample. In other 

10 words, the enhanced/manipulated spectral data more closely 
correlates with the actual spectrum emitted by the sample. 

The method comprises steps, some of which are not 
essential, which should be carried out on the digital 
spectral data, as follows. 

15 Step 1; To a raw spectrum (spectral or frequency- 

domain) defined by a set of 2 m pixels (where m is an integer 
number) , apply an Inverse Fast Fourier Transform (IFFT) . 2 m 
data points are required for applying Inverse Fast Fourier 
Transform algorithm. (A shifted IFFT provides. a result as a 

20 symmetrical pseudo-interf erogram; all computations are more 
elegant due to the symmetry) . An "interf erogram " type data 
set is obtained, in a 'pseudo'- time domain. The 
interferogram has M = 2 m (un-shifted IFFT) or 2 m +l (shifted 
IFFT) data points, depending on the way the dataset is 
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handled. Such an interf erogram 80 is shown in Figure 4. In 
this case, it has 2 m +l data points, it is symmetrical about 
time t=0 and has the time scale normalized to 1 seconds, 
from -0.5 to +0.5 seconds. 
5 Step 2 : Zero- fill the interf erogram thus adding 2 n -2 m 

data points with intensity equal to zero. This increases the 
number of data points to 2 n +l where n>m. The number 2 (n " m) 
gives the degree of Zero-filling. Zero-filling is a 
technique by which zero values are added (symmetrically only 

10 for the shifted IFFT) to the real and the imaginary part of 
the IFFT* for the new data points. In other words, the IFFT 
data between -0.5 and 0.5 on the pseudo-time scale remains 
unchanged and data with a value of zero is added to the IFFT 
between values of -1.0 to -0.5 and 0.5 to 1 (in this case 

15 where the time scale is enlarged by a factor of 2, the 
degree of zero-filling is 2) . 

Step 3: (optional, non-essential step) apply 
apodization to the interf erogram. Apodization is a 
multiplication of the imaginary and real part of the IFFT 

20 interferogram with a selected function in order to improve 
either the signal-to-noise ratio (equivalent to smoothing) , 
to the detriment of resolution, or to improve spectral 
resolution to the detriment of the signal-to-noise ratio. 
Examples of apodization are discussed below. 
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Step 4 ; apply a Fast Fourier Transform (FFT) to the 
results of steps 2 (or step 3, if used) . The obtained 
spectrum has 2 n points . 

As an example, Figures 5, 6, 7 and 8 show the effect of 
5 zero- filling to various degrees on an inverse Fourier 

Transform of a raw optical emission spectrum, followed by 
the FFT back to the wavelength domain. No apodization (i.e. 
step 3 outlined above has not been performed) has been 
performed on the data in the example shown in these figures. 
10 The raw spectrum was recorded with spectrograph 2 of ARL 
QUANTRIS using a CrNi steel sample to produce a spectrum. 

Figure 5 shows the raw data 90, obtained from the 
spectrometer, plotted as a function of pixel number against 
an arbitrary intensity value; the pixel-limited digital 
15 resolution is evident. Three examples of zero-filling are 
presented in figures 6, 7 and 8, where the factor of zero- 
filling used is 2, 4 and 8 respectively. 

As can be seen, zero-filling alone does not change 
spectral resolution (and indeed cannot do so since this is 
2 0 determined by the optical arrangement of the spectrometer) , 
but it does increase the digital resolution, and it also 
makes the shape of the peaks closer to the natural shape of 
the peaks (the manipulated data peaks having a "less- 
digitised look") . Minor details (previously hidden) are also 
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unveiled, such as the peak indicated by numeral 92 in 
figures 6, 7 and 8 which is not visible in the raw data. 
Furthermore, the raw data shows a spectral region 94 which 
comprises two adjacent spectral peaks 96, 98 of different 
5 intensity. However, the same region 100 in the 

manipulated/zero- filled data of figures 6, 7 or 8 shows two 

9 

peaks with substantially the same intensity. 

The resultant spectrum can be analysed to determine the 
peak position with an accuracy and precision which was 

10 previously not readily possible. Integration limits can be 
set with far greater precision (up to a factor 8 in the 
example shown in figure 8). Drift compensation can also be 
more precisely applied. This technique reduces the 
correction scale unit down to a logical pixel , instead of a 

15 physical pixel width. 

Comparison of data fitted using an embodiment of the present 
invention with known techniques. 

20 A comparison of data manipulated using an embodiment of 

the present invention with the technique of linear 
interpolation is now presented. Figure 9 shows some raw 
spectral data 100 presented as a histogram chart. The data 
comprises of a selection of twenty-six data points plotted 
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along the x-axis, each having different intensities which 
are represented by the height of each bar in the histogram, 
and which are plotted along the y-axis. It can be seen that 
the spectrum comprises two singlet peaks 102 and 104, and a 
5 doublet 106'. The first peak appears to have a FWHM value of 
between one and three pixel widths - it is certainly not 
possible to given an accurate estimation of the peak's FWHM 
value. Likewise with the other peaks shown in the spectrum. 
Referring to Figures 10 and 11, a plot of. the spectral 

10 data 100 of figure 9 is shown after the data has been 

subjected to linear interpolation of a fourth and eighth 
degree respectively. Briefly, a straight line is plotted 
between two' adjacent data points by plotting a number of 
contrived data points between real data points. If linear 

15. interpolation is carried out to two degrees, then two 

contrived data points are plotted between two adjacent real 
data points'. Likewise if linear interpolation is carried out 
the eighth degree, then eight contrived data points are 
plotted between two adjacent real data points. 

20 The results of linear interpolation plot 110 shown in 

figures 10 and 11 are plotted as histogram charts. It can be 
seen that linear interpolation only goes a small way to 
improving the raw data. For instance, it is still difficult 
to ascertain exactly the centre of the second peak 114, or 
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the centre of the second doublet peak 118. The situation is 
not improved by using a higher power linear interpolation 
plot 120 as shown in figure 11. By increasing the degree of 
interpolation it can be seen (from a comparison of figures 
5 10 and 11) that no great benefit is yielded for someone 

trying to determine the centre of peak 124, etc.: The centre 
of peak 124 appears to be in exactly the same position as 
the centre of peak 114. 

Referring now to figures 12 and 13 , a plot of the 

10 spectral data 100 of figure 9 is shown after the data has 

been subjected to manipulation according to an embodiment of 
the present invention of a fourth and eighth degree 
respectively without apodization. Once' again, figures 12 and 
13 show a selection of a spectrum after manipulating all 

15 8640 pixels in the spectrum, and not just the 26 pixels 
displayed. 

From the resultant manipulated data shown in figure 12 
and figure 13, it can be seen that it is much easier to 
discern details in the spectrum which were not otherwise 
20 apparent from either the raw data or the data which has 

undergone linear interpolation. For instance the centre of 
peak 134 and 144 can now be easily determined, compared to 
inadequate centre measurement of peaks 104, 114 or 124. 
Also, the shape of the doublet 13 6 and 14 6 is much better 
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defined with respect to the raw data and data having linear 
interpolation performed thereon. Features in the zero-filled 
data appear which were not apparent from either the raw data 
or the linearly interpolated data. For instance, the peak 
5 indicated by numerals 137 and 147 is not entirely apparent 
in either the raw data or data on which linear interpolation 
has been performed. 

Referring now to figure 14, a direct comparison of the 
plots shown in figures 11 and 13 is made. Here, the data 

10 from each of the respective data points is shown as a line 
plot, rather than a histogram. This makes the comparison 
easier to visualise. The data from figure 11 (that is, the 
linear interpolation to an eighth degree data) is indicated 
by line 150* (individual intensities represented by diamonds 

15 on line 150) . The data from figure 13 (that is the zero- 
filled data according to an embodiment of the present 
invention) is indicated by line 152 (individual intensities 
represented" by circles on line 152). 

The singlet peak 162 on the left hand side of the 

20 spectrum shows relatively good correlation between the zero 
filled data and the linear interpolation data. The FWHM of 
both sets of data are similar and both sets of data show 
good correlation for the predicted centre wavelength of this 
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peak. Also, the intensity of the peak 162 is similar for the 
linear interpolated and zero-filled data-- 

However, there are large differences in the 
characteristics of other peaks in the data which are readily 
5 apparent. For instance, the intensity of the peak 164 varies 
considerably between the two data sets. Likewise there is 
little correlation between the data sets with regards to the 
overall shape of the doublet 166, particularly the right- 
hand peak of the doublet. Also, the centre of the peaks is 

10 much clearer from the zero-filled data, particularly for the 
164 peak. It is not at all clear from the linear 
interpolated data exactly where the centre of peak 164 lies. 
It would appear to be anywhere between pixel number 19 to 20 
(on the x-axis) . However, the zero-filled data shows a clear 

15 discernable peak at one value. 

Apodization 

Apodization is a known method of further manipulating 
data to increase signal to noise ratio, to reduce artefacts, 
20 or to increase resolution. Essentially, apodization 

comprises imposing a function on the real and imaginary ■ data 
of a time-domain signal. Depending on the apodization 
function chosen, the resultant data can be further enhanced 
when it is transformed back into the wavelength domain. 
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Figure 15 shows various apodization functions: a so- 
called cosine square 180, a so-called shifted sine-bell 182 
and a Hamming function 184 . The cosine square function 
almost fits the signal 190 envelope which looks as a 
5 function centred at maximum, decaying monotonously towards 
the ends. The further the signal is from the centre of the 
graph, the more the intensity is reduced by multiplying by 
the function. The resultant signal envelope has a faster 
decay which corresponds in the spectral domain to a broader 

10 FWHM. Noise being constant over the interf erogram, the part • 
of the interf erogram with the worst signal to noise ratio is 
given less weight. In other words, the signal to noise ratio 
is improved at the expense of spectral resolution. Other 
functions having a monotonous decay (for instance an 

15 exponential function) perform in a similar way to the cosine 
square function. A variety of other known functions can be 
used to improve the resolution, for instance a shifted sine- 
bell or a Hamming function. Usually these functions put less 
weight on the part of the interf erogram 190 around time^Os, 

2 0 than on the intermediate parts, for example, around 
times=±0.10 - 0.30s. 

Figure 16, 17 and 18 show raw data, the effects of 
zero-filling only, and of apodization on a zero-filled 
spectrum, respectively. An Inverse Fast Fourier Transform is 
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applied to the raw optical emission spectrum 200 (Figure 
16) , followed by zero- filling (Figure 17) and by zero- 
filling and apodization (Figure 18) . Finally a Fast Fourier 
Transform returns the data to the wavelength domain. In this 
5 example a pure Al sample spectrum is used. Using as an 

apodization function cos(t) 2 (so-called square cosine bell, 
180 in figure 15) where t is the pseudo-time, a smoothing 
effect is obtained, as can be seen by comparing figures 17 
and 18. This smoothing is important for evaluating 

10 background regions and improving the signal-to-noise ratio 
for data at low concentration levels. It can be seen that 
the smoothing has improved the signal-to-noise ratio of the 
data set, but at the expense of spectral resolution; the 
line widths have increased. Other apodization functions are 

15 available that also improve the signal- tp-noise ratio in 
this way. Normally these are functions that make the 
envelope of the pseudo-time signal shorter. 

Improvements to the digital resolution are only 
worthwhile if the spectral resolution is greater than the 

20 digital resolution (since otherwise the former becomes the 
limiting factor) . Apodization techniques can be used to 
improve digital resolution or reduce artefacts in cases 
where the pseudo-time domain data is truncated (i.e. has a 
significant intensity at the extremities) . In such cases, 
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after Zero-filling, Fourier Transforming the signal leads to 
artefacts in the structure of the peaks or lines. When the 
intensity range of the spectrum (i.e., the dynamic range) is 
large, the artefacts arising from the largest peaks may have 
5 intensities comparable to those of the smallest peaks. To 
avoid this, .truncated signals in the pseudo-time domain can 
be apodized to more smoothly take the signal to zero. This 
apodization itself leads to line or peak broadening. To 
avoid this, linear prediction is applied to generate the 

10 additional data points required to . take the pseudo-time 
signal to zero, without introducing the line broadening • 
effects of the apodization function. 

Embodiments of the present invention thus improve 
spectral details and resolution, allow the use of cheaper. 

15 CCDs (possibly with fewer pixels than is presently required 
for necessary resolving powers) and/or to reduce the time 
taken to obtain accurate spectra. The time saving is 
particularly beneficial in sequential (scan) techniques. 

Although the invention has been described in connection 

20 with optical emission spectra it will be understood that the 
techniques are equally applicable to other forms of spectra. 
By way of example only, methods embodying the present 
invention can be applied to other spectrometers which 
produce intensity versus wavelength measurements, such as 
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Inductively Coupled Plasma-Optical Emission Spectroscopy, 
Energy Dispersive-X-Ray Fluorescence and Wavelength 
Dispersive-X-Ray Fluorescence. 

The foregoing discusses the application of an inverse 
5 Fourier Transform to data in the frequency or wavelength 
.domain, to produce a data set in what is referred to as the 
"pseudo-time domain". It is this pseudo-time domain data 
that is apodized and/or zero filled, before transformation 
back into the wavelength/frequency domain. It is however to 

10 be understood that the method described is. equally 

applicable to intensity data obtained as a function of mass 
(strictly, mass to charge ratio) . Such data is routinely 
obtained from, for example, Inductively Coupled Plasma-Mass 
Spectrometry, Gas Chromatography-Mass Spectrometry, organic 

15 MS-MS, Time of Flight (TOF) MS, or triple quadrupole 
techniques using, for example, electro-spray sources. 

In such cases, as with the intensity vs 
frequency/ wavelength embodiments described above, an inverse 
Fast Fourier Transform is first generated. The data thus 

2 0 transformed (in what is referred to here as the inverse 
transform mass domain) is apodized and zero filled as 
previously and then an FFT is applied to the resultant data 
to convert it back into the mass domain. 
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Figure 19a shows a "raw" mass spectrum in the absence 
of any further processing, on linear vertical and horizontal 
scales, and Figure 19b shows, the same "raw 11 mass spectrum on 
a log-linear scale. Figure 20a shows the same mass spectrum, 
5 again on a linear vertical and horizontal scale, but with 
the data post processed by applying an IFFT, Zero-filling to 
degree 4, and then applying an FFT to convert the data back 
into the mass domain. Finally Figure 20b shows the processed 
mass spectrum of Figure 20a but on a log linear scale. 

10 Two features are notable: firstly, the peaks in Figure 

2 0a and 2 0b are smoother than those in Figures 19a and 19b 
(exactly as with the earlier embodiments of Figures 16, 17 
and 18) , but secondly, a significant amount of artefacts are 
introduced away from the peaks, as is particularly apparent 

15 in the log linear plot of Figure 20b. This is a consequence 
of the nature of the mass spectrum in Figures 19a and 19b, 
which each represent an elemental mass spectrum consisting 
of a relatively small number of peaks across the measured 
mass range, with a large range of intensity values (a high 

2 0 dynamic range) . Artefacts due to the highest peaks are 

comparable to the smallest peaks of interest. Zero-filling 
in the regions that have no peaks leads to the artefacts 
seen best in Figure 2 0b. 
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By apodizing the data in the inverse transform mass 
domain prior to Zero- filling, a good deal of these artefacts 
can be removed. Figures 21a and 21b show the same processed 
mass spectrum as in Figures 20a and 20b, but this time with 
5 a cosine squared apodization function applied to the "raw" 
data from the mass spectrum once an IFFT has been applied to 
it to transform it into the inverse transform mass domain. 
The apodization is applied before the subsequent Zero- 
filling takes place, again in the inverse transform mass 

10 domain, and then an FFT back to the mass domain occurs. It 
will be seen by comparison of Figures 2 0b and 21b in 
particular that the number of artefacts away from the peaks 
have been reduced. Apodization using the cosine squared 
function has slightly broadened the peaks and slightly 

15 reduced the peak height but the noise has almost completely 
been removed. 

Apodization effectively acts to "weight" the Zero- 
filling preferentially around the peaks and so the choice of 
apodization function is dependent upon the spectral shape. 

2 0 Other functions might therefore be employed depending upon 
the anticipated nature of the mass spectrum. For example, 
whilst a cosine squared function appears to suit the 
concentrated peak profile of the elemental mass spectrum of 
Figures 19, 20 and 21, other functions might be even more 
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appropriate, and for other types of mass spectrum, still 
other apodization functions might be more suitable. For 
example, in so called MS -MS or MS n experiments carried out 
on organic molecules in triple quadrupole arrangements, a 
5 mass spectrum will often contain a continuous or almost 
continuous set of peaks over a selected mass range, as a 
result of the presence of multiple precursor and fragment 
ions. The technique described above is equally suitable to 
such spectra, and, due to the usual presence of a continuous 

10 • or quasi continuous set of peaks over a mass range with a 

relatively modest intensity (or dynamic) range, Zero-filling 
alone may produce a processed mass spectrum without 
significant noise added. 

In addition to ICP-MS and triple quadrupole 

15 applications, the skilled person will appreciate that the 
technique is equally applicable to mass spectra produced 
from still other mass spectrometer arrangements, including 
but not limited to magnetic sector devices, 3 -dimensional 
traps, time of flight (TOF) devices and the like. 

20 Moreover, imaging spectroscopes such as ICP-CID spectra 

(with a bi- dimensional FFT processing) can be treated by 
methods of this present invention. 

Furthermore, it can be applied to spectra that have 
been recorded by sequential spectrometers. In such cases, it 
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can save important scanning time (and costs) by increasing 
the scan step size by a factor of two or four, without 
prejudice to the final spectrum. 



